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Abstract 

In this paper, we construct a new class of blowup solutions with elementary functions to the 
3-dimensional compressible or incompressible Euler and Navier-Stokes equations. In detail, 
we obtain a class of global rotational exact solutions for the compressible fluids with 7 > 1: 

p = max 1 2_^L [c* 2 ^ x 2 + y 2 + z 2 - (xy + yz + xz)] - a(t)(x + y + z) + b(t)] , o| 7-1 

mi = a(t) + C (y - z) (1) 
u 2 = a(t) + C {-x + z) 
U3 = a(t) + C (x — y) 

where 

a(t) = co + at (2) 

and 

b(t) = 3c at + 2 C ?* 2 + C2 ( 3 ) 
with C, Co, ci and ci are arbitrary constants; 

And the corresponding blowup or global solutions for the incompressible Euler equations are 
also given. Our constructed solutions are similar to the famous Arnold-Beltrami-Childress 
(ABC) flow. The solutions with infinite energy can exhibit the interesting behaviors locally. 
Besides, the corresponding global solutions are also given for the compressible Euler equations. 
Furthermore, due to divu = for the solutions, the solutions also work for the 3-dimnsional 
incompressible Euler and Navier-Stokes equations. 
MSC: 35B40 35Q30 35C05 76D05 
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1 Introduction 

The iV-dimensional Euler equations can be formulated as the follows: 

p t +V ■ (pu) = 



p [u t + [u ■ V) tZ] + VP(p) = ^ 

where x — (x%, X2, xn) £ R , p = p(t, x) and u — u(t, x) — (ui,U2, ■■■u n ) £ R N are the density 
and the velocity respectively. The 7-law could be applied to the pressure function, i.e. 

P(p) = Kpi (5) 
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with K > and 7 > 1. The Euler equations @ can be rewritten by the scalar form, 



Ot t- — /ft:— ± C/Xfc ' ^ — 'ft: — J. uxfc ' 



i=l 



The Euler equations (jlj are the fundamental model in fluid mechanics [2] and [5]. 

Constructing exact solutions is a very important part in mathematical physics to understand 
the nonlinear behaviors of the system. For the pressureless fluids K — 0, a class of exact blowup 
solutions were given in Yuen and Yeung's papers [14] and |llj : 

/( W W'"--^ #2 > «. = gig (a; , + d . } f 0ri = l, 2 ,....,iV 

n adt) ,K ' (7) 

asj(*) = a i0 + ant 

with an arbitrary function / > and / € C 1 ; and cii, a^o > and an are constants. 
In particular, for a 2 < 0, the solutions blow up in the finite time T = —a\jai. There are other 
analytical solutions for the compressible Euler or Navier- Stokes equations in [9], jS], [6], [12], [13], 
[M] . [TT] . [15] and [16]. In addition, there are two papers to investigate for a special structural 
form which generalizes the exact solutions for Burgers' vortices in [7] and [3]. 

In this paper, we manipulate the elementary functions to construct some exact rotational 
solutions for the 3-dimensional compressible Euler equations (j4}: 

Theorem 1 For the 3-dimensional compressible Euler equations Q), there exists a class of rota- 
tional solutions: 
for 7 > 1 : 

p = max I Z=± [C 2 [x 2 +y 2 + z 2 - (xy + yz + xz)] - a(t)(x + y + z) + b(t)] , oX' 1 ' 1 

u 1 =a{t) + C(y-z) (8) 
U2 = a(t) + C (—x + z) 
u 3 = a(t) + C(x-y) 

where 

a(t) = c + cit (9) 

and 

3 

6(f) = 3c ci< + -eft 2 + c 2 (f 0) 

with C , Co, Ci and C2 are arbitrary constants; 
for 7 = 1; 

„C 2 [x 2 +y 2 +z 2 -(xy+yz+xz)]-a(t){x+y+z)+b(t) 

ui = a(t) + C (y — z) 
u 2 = a(t) + C (-x + z) 
u 3 = a(t) + C(x-y). 



(11) 



The solutions (0) and ill]) globally exist. 

Remark 2 In 1965, Arnold first introduced the famous Arnold- Beltrami- Childress (ABC) flow 

u\ = A sin z + C cos y 

U2 = B sinx + A cos z (12) 
U3 — C sin y + B cos x 



with constants A, B , C and a suitable pressure function P only for the incompressible Euler 
equations in JJJ. We observe that our solutions |3|) and {77]) are similar to the ABC flow. 
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Remark 3 The solutions with infinite energy can exhibit the interesting behaviors locally. The 
exact solutions with infinite energy of the systems may be regionally applicable to understand the 
great complexity that exists in turbulent phenomena. 

Remark 4 We notice that the rational functional form with a(t) =0 and C = 1 in solutions (0j 
and {Tip for the velocity u comes from Senba and Suzuki 's book MOtf . The velocities u in solutions 
|3p and 111]) are not spherically symmetric. 

Remark 5 The exact rotational solutions @] and ill]) could be good examples for testing numer- 
ical methods for fluid dynamics. 

2 Compressible Rotational Fluids 

The main technique of this article is just to use the primary assumption about the velocities u: 

ui = a(t) +C (y — z) 

u 2 = a(t) + C {-x + z) (13) 
"3 = a(t) + C (x - y) 

to substitute the governing equations (j4|) to construct the density function p to balance the system. 

Therefore, the proof is simple to be checked by direct computation: 

Proof. For the vacuum solutions p = 0, it is the trivial solutions for the system 

For 7 > 1, with non-vacuum solutions, we have the first momentum equation ©2,1: 

| Ul + Ul JL Ul + U2 JL Ul + U3 JL U1 + Kw-r-'JLp (i4) 

Odd d . . 

= -^l u l + u 2-^-Ui + U 3 — Ui -\ T^-P 1 (15) 

at ay oz 7 — 1 ox 

with 

-^-p 1 - 1 = C 2 (x 2 +y 2 + z 2 - (xy + yz + xz)) - a(t)(x + y + z) + b(t) (16) 
7 — 1 

which we could require the condition for the density function: 

/9 = max|^-i [C 2 [x 2 +y 2 + z 2 - (xy + yz + xz)] - a(t)(x + y + z) + b(t)] , oj^ , (17) 
to have 

= a(t) + [a(t) + C(-x + z)] — [o(t) + C{y-z)] + \a{t) +C(x-y)} — [o(t) + C(y- z)] 
d 

+ [C 2 [x 2 + y + z 2 - (xy + yz + xz)] - a{t)(x + y + z) + b(t)] (18) 
= a(t) + [a(t) + C{-x + z)\ C + [o(t) + C (x - y) (-£?)] 

+ [2C 2 x - C 2 y - C 2 z - a(t)] (19) 

= a(t) - 2C 2 x + C 2 z + C 2 y + 2C 2 x - C 2 y - C 2 z - a(t) (20) 

-0. (21) 

Similarly by the nice symmetry of the functions ([8jl, we could balance the second momentum 
equation ©2,2: 

d d d d 1 , 

■7T7W2 +U1—U2 +U 3 —U 2 H T^rP ( 22 ) 

at ox oz 7 - 1 oy 
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a(t) + [a(t) + C(y-z)} — [a(t) + C (-x + z)] + [a(t) + C(x-y)} — [a(t) + C (-x + z)] 



d 

+ ^ [C 2 [(x 2 + y 2 + z 2 ) - (xy + yz + xz)} - a(t)(x + y + z) + 6(f)] 

= a(t) + [a(t) +C(y- z)} (-C) + [o(t) + C(x-y)]C 
+ 2C 2 y - C 2 x - C 2 z - a(t) 

= 0. 



(23) 
(24) 

(25) 



For the last equation d4} 2 ,3, we have 





d d K^f d 

u 3 + ui—u 3 + u 2 tt u 3 H T^rP 

at ox ay 7 — 1 az 



(26) 



= d(t) + [a(f) + C(y-z)] — [a(t) +C(x-y)} + [a(t) + C {-x + z)] — [a{t) + C(x- y)} 

Q 

+ g- [° 2 [( x +y 2 + z2 ) - ( x y + y z + xz )\ - ^)( x + y + z ) + &(*)] ( 27 ) 

= d(t) + [a(t) + C(y-z)]C+ [a(t) + C {-x + z)] (-C) + 2C 2 z - C 2 y - C 2 x - a(t) (28) 

= 0. (29) 
For the mass equation ©1, we have with the density (fl~7f to obtain: 

p t + V • Con) (30) 



+ 



at 

<9x 
c)y 

dz 



1 — 1 

(C 2 (a; 2 + y 2 + z 2 - [xy + yz + xz)) - a(t)(x + y + z) + b(t)) 

7 ~ ^ ^2 1 „,2 1 _2 



(C 2 (x 2 +y 2 + z 2 - (xy + yz + xz)) - a(t)(x + y + z) + b(t)) 



T — 1 

(C 2 (x 2 + y 2 + z 2 - (xy + yz + xz)) - d(i)(x + y + z) + b(t)) 
(C 2 (x 2 + y 2 + z 2 - (xy + yz + xz)) - a(t)(x + y + z) + b(t)) 



(a(t) + C(y - z)) 

(a(t) + C (-x + z)) 

(a(t) + C (x - y)) 
(31) 



7-1 



7 — 1 

(C 2 (x 2 + y 2 + z 2 - (xy + yz + xz)) - a(t)(x + y + z) + b(t)) 



7-1 
if 7 



-S(<)(x + y + z) + 6(f) 
+ [(2C 2 x - C 2 y - C 2 z - a(t)) (a(t) + C(y - z))] 

[(2C 2 y - C 2 x - C 2 z - a(t)) (a(t) + C (-x + z))] 
+ [(2C 2 z - C 2 x - C 2 y - a(t)) (a(t) + C(x - y))] 



(32) 



1 



1 — 1 

(C 2 (x 2 +y 2 + z 2 - (xy + yz + xz)) - a(t)(x +y + z) + b(t)) 



-a(t)(x + y + z) + 6(t) 
+ (2C 2 x-C 2 y-C 2 z-d(t))a(t)+C(2C 2 x-C 2 y-C 2 z-d(t))(y-z) , 
+ (2C 2 y - C 2 x - C 2 z - a(t)) a(t) + C (2C 2 y - C 2 x - C 2 z - a(t)) (-x + z) ' ' 
+ (2C 2 z - C 2 x - C* 2 y - d(t)) a(<) + C (2C 2 z - C 2 x - C* 2 y - a(t)) (x - y) 
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(C 2 (x 2 +y 2 + z 2 - (xy + yz + xz)) - a(t)(x + y + z) + b(t)) 



j-a(i)(a; +y + z) + b(t) - 3d(t)a(i)} 



by choosing 

a{t) — c + cit 

and 

b(t) = 3a(t)a(t) 

b(t) = 3ci (co + cit) 
3 

b(t) = 3c oCl t + -c\t 2 + c 2 

where c 2 is the arbitrary constant. 

As for 7=1, the proof is similar, we may skip the details here. 
It is clear to see that the solutions © and (fTTj) globally exist. 
The proof is completed. ■ 

Here the masses of the solutions ((8} and (fTTj) are infinite: 



pdx = +oo. 



(34) 
(35) 

(36) 

(37) 
(38) 

(39) 



(40) 



3 Incompressible Rotational Fluids 

For the 3-dimensional incompressible Euler equations: 



divu = 
p [ut + (u ■ V) u] + VP(p) = 0, 



(41) 



the solutions with elliptical symmetric velocity in [5], [5] and [TB] for the compressible Euler 
equations could not be applied to the incompressible Euler equations (|4"Tj) because of the 
linear functional form: 

ui = ci(t)x + di(t) 

u 2 = c 2 {t)y + d 2 (t) (42) 
u-i = c 3 (t)z + d 3 (t) 



with some functions c,(i) and di(t) to have 



div u 7^ 



(43) 



for non-trivial solutions c\(t) + c 2 {t) + Ca(t) ^ 0. However, our solutions ([8]) and ([TT]) could be 
applied to the incompressible Euler equations (14T1) . The set of the corresponding solutions for the 
incompressible Euler equations (|4"Tj) could be larger than the compressible ones (HJ): 

Theorem 6 For the 3-dimensional incompressible Euler equations there exists a class of 

rotational solutions: 



P(p) 
p 



C 2 [x 2 + y 2 + z 2 - (xy + yz + xz)] - [oi(t) + C(a 2 {t) + a 3 (t))} x 
- [aa(t) + C(oi(i) + a 3 (t))] y - [d a (t) + C(ai(t) + a 2 {t))] z + b(t) 

ui = ai(t) + C(y - z) 
u 2 = a 2 (t) + C ( — x + z) 

u-3 = a 3 (t) +C(x-y) 



(44) 



() 
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P(p) = 
p 



and 

-C' 2 (x 2 + y 2 + z 2 + xy + yz + xz) 
-[a 1 {t) + C(a 2 (t) + a 3 {t))]x 
- [a2(t) + C(ai(t) + a 3 (t))]y 
- [d 3 (t) + C(ai(t) + 02(f))] z + b(t) J (45) 
ui = a x (t)+C(y + z) 
u 2 = a 2 (t) + C (x + z) 
u 3 = 03(f) + C (x + y) 

where Oj(f) is an arbitrary local C 1 function for i — 1, 2 or 3, b(t) is an arbitrary function, C is 
an arbitrary constant. 
In particular, 

1. if |a;(T)| = 00 or \cii(T)\ — 00 with the first finite positive constant T , the solutions {44}) and 
J^5| ) blow up in the finite time T ; 

2. For global C 1 functions ai(t), the solutions and J^5[ ) globally exist. 

Proof. For the checking of the functions (|4"4")l . the detail is similar to the proof of Theorem 1: 
For the checking of the functions (|45|) . we have for the first momentum equation (|4Tj) 9 9: 



d d d d d (P 

—Ui + Ui—Ul + + U3—UX + — — 

at ox ay az ox \p 



d d d d (P 

TT7 U 1 + U 2^-"l + U 3 — Ul + — — 

at ay az ox \ p 



by requiring the pressure function 

P 
P 

to have 



C 2 [x 2 +y 2 + z 2 - (xy + yz + xz)] - [di(t) + C(a 2 (t) + a 3 (t))] x 
- [a 2 (t) + C( ai (t) + oa(t))] y - [6 s (t) + C( ai (t) + a 2 (t))] z + b(t) 



(46) 
(47) 

(48) 



= 01 (t) + [a 2 (t) + C(-x + z)] — [a x (t) +C(y-z)] + [a 3 (t) +C(x-y)\ — [a x (t) + C (y - z)] 

ay az 



d_ 

dx 



C 2 [x 2 +y 2 + z 2 - (xy + yz + xz)] - [ai(t) + C(a 2 (t) + a 3 (t))] x 
- [6 a (t) + C(a x (t) + a 3 (t))] y - [03(f) + C(oi(t) + a 2 (f))] z + 6(f) 



(49) 



= ai(f) + [02 (t) + C(-a: + z)]C + [a 3 (t) +C(x~y) (-C)] 

+ [2C 2 x - C 2 y - C 2 z - d(f) - C(o 2 (f) + o 3 (f))] (50) 

= di(f) + C(a 2 (t) + a 3 (f)) - 2C 2 x + C 2 z + C 2 y + 2C 2 x - C 2 y - C 2 z - d(f) - C(a 2 (t) + a 3 (t)) (51) 

= 0. (52) 



Similarly by the nice symmetry of the functions (|44[) . we could balance the second momentum 
equation (gT^ 2 : 

d d 8 d (P\ 

-ttU 2 + ui—u 2 + u 3 —u 2 + — — (53) 
at Ox az ay \ p J 



s-1 Pi 

a 2 (t) + [01 (t) +C(y-z)] — [o 2 (f) + C + «)] + [03(f) + C (x - y)] — [a 2 (t) + C (-x + z)] 



d_ 

dy 



C 2 [x 2 +y 2 + z 2 - (xy + yz + xz)] - [di(f) + C(a 2 (t) + a 3 (t))] x 
- [a 2 (t) + C(oi(f) + o 3 (f))] y - [d 3 (f) + C(oi(f) + o 2 (f))] z + b(t) 



(54) 



= a 2 (t) + [01 (t) + C (y - z)] (-C) + [03(f) + C(i - y)] C 
2C 2 y - C 2 x - C 2 z - a 2 (t) - C(ai(f) + 03(f)) 



(55) 
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0. 



(56) 



For the last equation pljb q. we have 



d (P 



d d d 

+ Ut—Ua + U 2 —U 3 + — I — 

at ox oy Oz \ p 



(57) 



= 63 (t) + [01(f) + C (y - z)] — [a 3 (t) +C(x-y)} + [a 2 (t) + C (-a; + z)] — [a 3 (t) + C (x - y)] 



+ 



d_ 

dz 



C 2 [x 2 +y 2 + z 2 - (xy + yz + xz)] - [d x (t) + C(a 2 (t) + a 3 (t))] x 
- [a 2 (t) + C(oi(t) + a 3 (t))] y - [03 (t) + C(ai(i) + 03(f))] 2 + &(t) 



(58) 



= d 3 (f)+[a x (f) + C (y - z)] C+ [02(f) + C (-x + z)\ {-C)+2C 2 z~C 2 y-C 2 x-a 3 {t)-C{ ai (t)+a 2 {t)) 

(59) 

= 0. (60) 
Next, for the checking of the solutions (|45p . we have for the first momentum equation (|4"Tj) 9 1 : 



d 



d 



d 



d (P 



—ui + u 2 —ui + u 3 tt u i + s~ \ — 
at ay oz ox \ p 



(61) 



with the pressure function 

P 
P 

to have 



-C 2 (x 2 + y 2 + z 2 + xy + yz + xz) - [di(f) + C(a 2 (t) + a 3 (t))\ x 
- [d 2 (f) + C(oi(t) + o 3 (f))] y - [d 3 (f) + C(ai(t) + o a (t))] « + 6(t). 



(62) 



= di (f) + [o 2 (f) + C(x + z)] — [ 0l (t) + C (y + z)] + [03 (t) + C (a; + y)] ^ [<n(t) + C (y + z)] 



_9_ 

dx 



-C 2 (x 2 + y 2 + z 2 + xy + yz + xz) - [di(t) + C(a 2 (f) + 03(f))] x 
- [d 2 (f) + C(oi(t) + 03 (t))] y - [d s (t) + C(ai(t) + o 2 (f))] z + &(f) 



(63) 



= 01 (t) + [02(f) + C(s + z)] C+ [03(f) + C (x + y)] C - 2C 2 x - C 2 y - C 2 z - d(f) - C(o 2 (t) + a 3 (f)) 

(64) 

= di (i) + C(a 2 (t) + a 3 (f)) + 2C 2 x + C 2 z + C 2 y - 2C 2 x - C 2 y - C* 2 z - di (f) - C*(a 2 (t) + a 8 (t)) (65) 

= 0. (66) 



Similarly, we can balance the second momentum equation (|41|>9 9 with the symmetric form of the 
functions (1431) : 

d d d d [P\ , . 

= d-t U2 + Ul d-x U2 + U3 d-z U2+ dy-{^) (67) 

d d d f P 
= d 2 (f)+[oi(f) +C{y + z)] — [o 2 (f) + C (x + z)] + [a 3 (<) + C (x + y)] — [a 2 {t) + C (x + z)] + — I - 

ox oz oy \ p 



a 2 {t) + [a 1 {t) + C(y + z)]C+[a 3 {t) + C{x + y)]C-2C 2 y~C 2 x-C 2 z-a 2 (t)- C(oi (i) + a 3 (f )) 

= 0. (69) 



For the last equation plj) ? 3, we could get 



d 



d 



d (P 



-U 3 + U1—U3 + U 2 —U 3 + „ , 

ot ox oy oz \ p 



(70) 



o 3 (t)+[oi(t) + C (y + z)] — [a 3 (t) + C (a: + y)]+[aa(t) + C (1 + z)] — [a 3 (t) + C (x + y)}+— ( - 
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= a 3 (t) + [ai (t) + C {y + z)] C + [a 2 (t) + C (x + z)] C - 2C 2 z -C 2 y-C 2 x-a 3 (i) - C(ai (f ) + a 2 (*) ) 

(71) 

= 0. (72) 

It is clear to see that 

1. if |<Zi(T)| = oo or |d,-(T)| = oo for i = 1, 2 or 3, with the first finite positive constant T, the 
solutions (|44|) and (|45|) blow up in the finite time T; 

2. for all global C 1 functions aj(i), the solutions (|44|) and (|45|) globally exist. 
The proof is completed. ■ 

Here the kinetic energy of the solutions (|44|) and (j45"T) for the incompressible fluids are infinite: 

1 

2 



u = +oo. (73) 

And the blowup solutions (|44l) and (1451) are the other examples with infinite energy for the incom- 
pressible Euler equations (1411) . with respect to the blowup cylindrical ones in Gibbon, Moore and 
Stuart's work [3]. 

Additionally it is because 

Au = (74) 

in the solutions ([HJ and (fTT|) to be the corresponding solutions for the 3-dimensional compressible 
Navier-Stokes equations: 

p t +V ■ (pu) = , . 

p[u t + (u-V)u] + yP(p) = nAu K } 

with \i > 0. 

And the solutions (|4"4"|) and (|43|) are the corresponding ones for the 3-dimensional incompressible 
Navier-Stokes equations: 

divu = , , 

(76) 



p [u t + (u-V)u\ + VP(p) = /xAu 



with /i > 0. 



4 Discussion 

1. We observe that the constructed density functions p (© and (fTTj) . and the pressure functions 
(HU) and (|43|) share some nice algebraic geometry, but we do not understand for the implication 
of this geometric structure for the fluids. Could the solutions with small perturbations of these 
solutions converge to the original ones? 

2. Could we just cut a finite partial volume of the solutions to construct new week solutions? 
If we could, what happens for the corresponding weak solutions of the classical ones with infinite 
energy? 

3. By taking consideration with Makino's blowup solutions form (|42l) . we could guess that the 
more general functional form for the 3-dimensional velocities u in Cartesian coordinate: 

ui = ai(t) + bu(t)x + b 12 (t)y + b lz (t)z 

« 2 = em (t) + b 21 (t)x + b 22 {i)y + b 23 {t)z (77) 
u 3 = a 3 (t) + b 31 (t)x + b 32 {t)y + b 33 (t)z 

with the local smooth time C 1 functions a,i(t) and b{j(t) for i, j = 1, 2 and 3, 
could be the more general solutions structure for the systems. Then, it is more difficult to handle 
the linear functional structure for the higher dimensional cases R N with N > 4. Could we construct 
a blowup example with rotation for the incompressible fluids? In principle, we could adopt the 
guessing approach to construct the density functions. Some trials with luck are needed to obtain 
a right novel functional form to construct solutions for the systems, but it costs a lot of time 
for computation. Thus, we hope to arouse the readers' interests to investigate more systematic 
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approaches to calculate the possible particular solutions with the linear functional velocities u (|77p 
in Cartesian coordinate or other coordinates. 

Finally, further researches are needed for understanding this class ([77)) of blowup or global 
solutions for the systems in the future. 
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